Two-fluid hydrodynamics for a trapped weakly-interacting Bose gas 
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We derive the coupled equations of motion for the condensate (superfluid) and non-condensate 
(normal fluid) degrees of freedom in a trapped Bose gas at finite temperatures. Our results are 
based on the Hartree-Fock-Popov approximation for the time-dependent condensate wavefunction 
and an assumption of local equilibrium for the non-condensate atoms. In the case of a uniform 
weakly-interacting gas, our formalism gives a microscopic derivation of the well-known two-fluid 
equations of Landau. The collective modes in a parabolically trapped Bose gas include the analogue 
of the out-of-phase second sound mode in uniform systems. 
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The low frequency dynamics of superfluid *He is com- 
monly described using the two-fluid phenomenology first 
developed by Tisza and Landau This descrip- 
tion, later shown to be a consequence of a Bose broken 
symmetry ||,^, accounts for the characteristic features 
associated with superfluidity in terms of the relative mo- 
tion of normal fluid and superfluid degrees of freedom. 
In particular, it predicts the existence of second sound 
as an out-of-phase oscillation of the two components. 
In the present letter, we give a microscopic derivation 
of the analogous two-fluid equations for an inhomoge- 
neous weakly-interacting gas of trapped atoms. In this 
situation, the superfluid is identified with the conden- 
sate atoms as described by a macroscopic wavefunction, 
while the normal fluid density is associated with the non- 
condensate thermal cloud. In the uniform density limit, 
we show that our equations are consistent with the stan- 
dard two-fluid equations H,^. However, the hydrody- 
namic behavior of Bose gases is quite different from that 
of a Bose liquid such as superfluid ^He. 

Our analysis is based on the equation of motion of 
the macroscopic Bose wavefunction $(r, i), as de- 

termined within the time-dependent Hartree-Fock-Popov 
(HFP) approximation, which is a generalization of Refs. 

and 1^. This condensate wavefunction is coupled 
to the thermally excited atoms making up the non- 
condensate which we assume is described by a semi- 
classical phase space distribution function /(r, p, t). The 
further assumption that collisions are sufficiently rapid 
to force local equilibrium within the thermal cloud leads 
to a set of hydrodynamic conservation laws for the non- 
condensate component. Our final result is a set of equa- 
tions for the condensate and non-condensate degrees of 
freedom which extends recent work on the hydrodynam- 
ics of trapped Bose gases to below Tbec- 



We first consider the dynamics of the condensate. 
As usual, the Bose field operator is conveniently sepa- 
rated into condensate and non-condensate parts: "0(1") = 
$(r) -|- V'(r)- For an arbitrary nonequilibrium state, 
the spatially and time- varying macroscopic wavefunction 
$(r,t) = {xp{T))t is described within the time-dependent 
Hartree-Fock-Popov approximation by the equation of 
motion 



5$(r, t) 
dt 



2m 



+ Uextir) + 2gh{r, t) + gndr, t) 



X $(r,t) = ■H{v,t)^{v,t) . 



(1) 



Here the nonequilibrium non-condensate density is given 
by h{r,t) = {^jj'' {r)i(;{r))t, the condensate density is 
nc{r,t) — |$(r, t)p and g = ATrah^/m is the interac- 
tion strength. Eq.(|l|) represents a natural extension of 
recent work which approximated n(r, t) in (|^) by 
the equilibrium value no(r), thereby ignoring the collec- 
tive behavior of the non-condensate. In contrast, the 
time-dependent condensate wavefunction <i>(r,t) in (^ is 
coupled into the fluctuations of the non-condensate and 
a dynamical equation for the latter is also required. 

It is convenient to recast the condensate equation 
of motion into a pair of hydrodynamic equations us- 
ing the amplitude and phase representation $(r, i) = 
'I>(r,t)|e'^('"'*l Substituting this form into (Q) and sepa- 
rating real and imaginary parts, one finds 

f \ 
dt \ c cj 



dt 



1. 



= -Va 



(2) 



where Vc(r, t) 
defined by 



fi'Vd{r,t)/m. The potential 0(r,t) is 



1 



Mr,t)\ 



H(r,i)|$(r,t)| 



(3) 



where 7l!(r,t) is the HFP Hamiltonian given in (|^). An- 
ticipating the identification of Vc with the superfluid ve- 
locity V5, we see that ip plays the role of the chemical 
potential associated with the superfluid motion ||^. 

We next consider the dynamics of the non-condensate 
in the low-frequency collision-dominated hydrodynamic 
regime. In the semi-classical limit valid at finite tem- 
peratures 1^ (with fcsT ^ hujo, gno{r), where loq is a 
characteristic trap frequency), the dynamics can be for- 
mulated in terms of a quantum kinetic equation for the 
distribution function /(r, p, i) |^ 10|. Since (|^) implies 
that the number of particles in the condensate is con- 
served, we must for consistency exclude those processes 
which scatter atoms in and out of the condensate . In 
this situation, only collisions between excited atoms are 
relevant and wc can use the kinetic equation 



d_ 

di 



Vr-VU{r,t)-V, 



/(r,P,i) = 



a/ 

dt 



(4) 



Here U{r,t) = Uext{j^) + '^g['n'{^,t) + nc{r,t)] includes the 
self-consistent Hartree-Fock dynamic mean field in which 
the condensate part 2gnc{r,t) can be viewed as an addi- 
tional external field acting on the non-condensate. 

The required hydrodynamic equations are obtained 
from by making the further assumption that colli- 
sions enforce the distribution function to take the local 
equilibrium form M 



/o(r,p,t) 



1 



exp [f3{^[p - mv„]2 + U - n}] - I 



(5) 



where the thermodynamic variables /3, v„ and fi, together 
with U, all depend on r and t. When (||) is substituted 
into (^, the collision integral vanishes by virtue of the 
local equilibrium form of the distribution function. Tak- 
ing moments of the resulting equation with respect to 1, 
Pfj_ and , we obtain the closed set of equations |jll|] 



9n „ /- 
— + V • nv„ 

at 







+ (v„ • V)v„ 

at 



de 
di 



= - VP- fiVU 
V • (ev„) = v„ • VP . 



The non-condensate density n(r, t) is given by 

f d^p 1 
n{r,t)= J -^/o(r,p,t) = -^5r3/2(^(r,i)) 



(6) 



(7) 



with z{r,t) = e/'("-,t)[M(r,t)-c/(r,i)] ^^(r,^) = {27rh'^ / 
TOfcsr(r, t))^/2. The quantity e is the non-convective 
part of the kinetic energy density defined with v,i — 



m 



^ Similarly, P = |e is the kinetic contribution to 
the local equilibrium pressure defined by 



1 



.95/2(^(1-, t)). 



(8) 



Eqs. (g) and (^) constitute our full set of nonlinear hy- 
drodynamic equations for a trapped Bose gas at finite 
temperatures. 

At the level of approximation we are considering, the 
condensate and non-condensate satisfy separate continu- 
ity equations. Combining these two equations gives the 
expected two-fiuid continuity equation 



dn 
'dt 



(9) 



where n = h + ric and j = nv„ -|- ricVc- We identify the 
normal fluid density with n(r, t) and the superfluid den- 
sity with ric(r, t). The former identification is supported 
by noting that (|^) can be expressed equivalently as 



n(r, t) 



d^p p^ dfo{ep) 



(10) 



where ep{r,t) = ^ + U{r,t) — /i(r,t) is the excitation 
energy. This is the usual Landau formula for the normal 
fluid density §,|. 

The linearized version of (||) and (^ allows one to 
consider small amplitude oscillations about equilibrium. 
The equilibrium condensate wavefunction is determined 
by the solution of 



Ho(r)$o(r) 



Uext{r) + 2gho{r) 



2m 



gncoir) $o(r) = Mo$o(i-) 



(11) 



with nco(r) — \^o{r)\'^. The equilibrium non-condensate 
density no{r) is given by ^ with the equilibrium fugacity 
deflned as za = e''o['^°-'^o('')l , with Uo{r) = ?7ext(r) -t- 
2gno{r). Eq.(|ll|) and the equihbrium version of (0) must 
be solved self-consistently. 

The linearization of (^ around equilibrium leads to 
the condensate equations 



dSric 
dSvc 



-VScj), 



(12) 



where 



l'i>o(r)| 



[Ho{r) - ^io]Smr,t)\ 

+ gSndr, t) + 2g5n{v, t) . (13) 
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In arriving at this result, we have noted that V(/)o(r) = 
in equilibrium. Similarly, the linearization of (^) leads to 
the equations 



d5fi 
'di 

dt 

ddP 



= -VSP - SnVUo - 2gnoV{Sh 

5 2 
= -- V • (Po^v„) + -(5v„ • VPo , 



(14) 



where Pq (r) is the equilibrium kinetic pressure which sat- 
isfies VPo = — noVf/o- Above Tbec, these equations 
reduce to those of Ref. [Q if we ignore the effect of inter- 
actions {g = 0). 

For its intrinsic interest, and in order to better under- 
stand the implications of our two-fluid equations, we now 
consider the limit of a homogeneous system {Uextii^) ~ 0). 
In this special case, (11 ) yields a uniform condensate with 
the chemical potential having the Thomas-Fermi (TF) 
value /io — 2f/no + gnco- Taking all equilibrium quan- 
tities to be spatially-independent, and noting that the 
first term on the right hand side of (|l3|) can be neglected 
in the long- wavelength limit, the two velocity equations 
reduce to 



m- 



dSvc 
dt 

- dSVn 

mnQ- 



dt 



with 



-2gV(5n — g'VSric 

-VSP ~ 2gnoVSn - 2gfio'Vdnc 



dSP ^ 
^^-gPoV.Jv,, 



(15) 



(16) 



The two equations in (|1^) can be combined to give 
25("c0 + no)V5n - g{nco 



-V6P 
-VSP 



2no) V^Uc 
(17) 



One can verify that ( p7| ) is consistent with the following 
expression for the total local thermodynamic pressure 



P = P 



2nh 



(18) 



which is the equation of state at the level of approxima- 
tion we are considering. Within the same approximation, 
the total internal energy density is given by 



e(r,t) = e{r,t) + ig(V-t(r)V-t(r)?A(r)^(r))t 
~ e(r, t) -\- ^g[n'^ + "^nn - n^] . 



(19) 



These results and the equilibrium thermodynamic re- 
lation 



P = sT 



fin 



(20) 



allow us to identify the entropy density s. Using /io 
g(2no + n^o) together with (fsl) and (|l|), (M) gives 



5 ~ 



gricona . 



(21) 



Here Po is the equilibrium kinetic pressure defined in (jg) , 
with z = zo = e/SoCpn-Sffno) ^ g-/3gn,o. Using the local 
equilibrium expression for P in (^), one finds that the 
fluctuation in the total pressure P is given by 



SP — sqST + fio^/i + gnco{2Sn -\- Suc) , 



(22) 



with So defined by (|21|). Comparing this to the thermo- 
dynamic relation SP — sqST + ngSfi, we arrive at 



Sfi = S{2gn -\- guc) 



(23) 



This result confirms (in the case of a uniform equilibrium 
density) that S(j){r, t) in (|l2|) and ( p^ ) is indeed the fluc- 
tuation in the local chemical potential. More generally, 
we have verifled that (^) is equivalent to the key Landau 
equation for superfluid flow 



dt 



= -v^i{v,t) . 



(24) 



Finally, using (y_9|) we flnd 
3 - 

Se — + .9'^co(2(^n- + Sric) -f 2ghQSn . (25) 

Inserting this result into the thermodynamic relation 
Se = TqSs -t- fioSn^ a simple calculation gives 



ToSs = -SP + gncoSn 



(26) 



Taking the time-derivative, and using (|l^ ) and the conti- 
nuity equation for Sn(r,t), we obtain the linearized form 
of Landau's entropy conservation equation 



dSs 
'dt 



-soV ■ (5v„ 



(27) 



In the uniform case, Eqs. ( [l5| ) and (^6|) are easily 
solved to give the expected flrst and second sound phonon 
modes. In contrast to superfluid "^He, second sound in 
a gas involves a condensate oscillation largely uncou- 
pled from the non-condensate, with a velocity given by 
U2 = [grico/mY/'^. 

As a speciflc application to a trapped Bose gas, we con- 
sider the center-of-mass mode solution for an anisotropic 
parabohc potential Uext{i^) = \^{^x^'^ + ^lu"^ +^1^^)- 
On the basis of the generalized Kohn theorem |]l2|, one 
expects three modes in which the gas oscillates rigidly 
along each of the principal directions at the appropri- 
ate trap frequency cjj . Denoting the displacement of the 
gas by r]{t), the condensate and non-condensate densi- 
ties both behave as no(r — vii))^ giving a density fluc- 
tuation Sn{r,t) = — Vno(r) • ri{t), and a velocity field 



3 



v(r, = ri{t) which is spatially independent. One can 
easily verify that our linearized equations (|lj) and (|lj) 
admit a solution of this kind with both the conden- 
sate and non-condensate having identical displacements 
T](t) = rjoCOSLUit. Thus, in contrast to the static HFP 
used in Ref . |5| , the dynamic HFP theory given by (|^) is 
consistent with the generalized Kohn theorem. 

An estimate of the frequency of the out-of-phase dipole 
mode in a trapped gas can be based on the intuitive 
idea that the condensate and non-condensate oscillate 
rigidly against each other. The problem is then equiva- 
lent to a two-particle problem with masses Mc — raNc 
and Mn = mN , each confined in the parabolic poten- 
tial and coupled together by a spring with force constant 
k. Since the interaction energy between the two compo- 
nents is Eint =2(7/ (irnc(r)n(r), k can be determined by 
considering small displacements of the two components 
along the i-th direction. The equations of motion for the 
two coupled masses can then be solved, giving the in- 
phase center-of-mass mode at the frequency Ui discussed 
above, and an out-of-phase mode at the frequency 



25 



Mc 4- M„ 

MrM„. 



dr 



ducoir) 9no(r) 



(28) 



This result is confirmed by a more detailed treatment 
p3[ in which the hydrodynamic equations (p^-p^) are 
cast into the form of a variational principle. In Fig. 1 
we show the in- and out-of-phase mode frequencies as a 
function of temperature for a gas trapped in an isotropic 
parabolic potential. We also show the relative ampli- 
tudes of the condensate and non-condensate oscillations 
which satisfy the condition Mcfjc + Af„r7„ = 0, corre- 
sponding to the center of mass being stationary for any 
temperature. A measurement of these amplitudes would 
therefore determine directly the ratio Nc/N of the two 
fluid components. This mode of the trapped Bose gas is 
the analogue of the usual second sound mode ^ in bulk 
superfluid ^He (for which psvg -I- pn'Vn — 0). 

In summary, we have given a microscopic derivation 
of the hydrodynamic two-fluid equations for a trapped 
weakly-interacting Bose g summarized by Eqs. ( |l^ ) 
- ( p^ ) . Our analysis shows how such a two-fluid descrip- 
tion arises naturally from the existence of a macroscopic 
condensate wavefunction which is coupled to the non- 
condensate atoms. We have also shown explicitly how our 
linearized two-fluid equations reduce to the well-known 
equations (|), (|), and ^) of Landau [|§ when 
applied to a uniform Bosc-condcnsed gas. Finally, we 
have shown that our equations are consistent with the 
generalized Kohn theorem and obtained the out-of-phase 
mode analogue of second-sound in a trapped Bose gas. A 
more complete discussion of the two-fluid dynamics will 
be given elsewhere . 
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FIG. 1. (a) Mode frequencies for the in-phase (solid dots) 
and out-of-phase (open dots) dipole modes vs. temperature 
for 2000 Rb atoms in an isotropic parabolic trap (see Ref. [5] 
for values of the physical parameters used), (b) Condensate 
(r^c) and non-condensate (r;„) amplitudes for the out-of-phase 
dipole mode, (c) Fraction of atoms in the condensate as a 
function of temperature. 
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